Resistive axisymmetric equilibria with arbitrary flow 
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An analysis of axisymmetric equilibria with arbitrary incompressible flow and finite resistivity 
is presented. It is shown that with large aspect ratio approximation or vanishing poloidal cur- 
rent, a uniform conductivity profile is consistent with equilibrium flows. Also a comment made on 
coexistence of both toroidal and poloidal flows in an axisymmetric field-reversed configuration. 
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Calculating the equilibrium is one of the fundamental 
problems of magnetically confined plasmas. Most stud- 
ies are directed toward finding an ideal (i.e. infinitely 
conducting) magnetohydrodynamic (MHD) equilibria in 
an axisymmetxic plasma. The- earliest calculations are 
those of GracEl and Shafrajpovo, leading to the famous 
Grad-Shafranov equatioinu. The ideal and static Grad- 
Shafranov equation is an elliptic differential equation in 
the magnetic flux function ip with two arbitrary surface 
quantities as the pressure p(ip) and the poloidal current 

Consequently, there have been attempts to include var- 
ious effects e.g. mass flow into the equilibrium equations. 
An equivalent of Grad-Shafranov equation in an axisym- 
metric ideaL plasma with arbitrary flow has been, given 
by Hameiricl. In some recent works, SteinhauerQ deals 
with a generalization of Grad-Shafranov equilibria in p. 
multi-fluid with flow and Throumoulopoulos and TassoEl 
consider a helically symmetric equilibria with flow. The 
situation with flow becomes more realistic when one re- 
alizes the existence of equilibrium flows both in toroidal 
and poloidal directions in tokamaks following momentum 
deposition through heating by neutral beam injectionLTu. 
With equilibrium flows, the resultant governing differen- 
tial equation does not remain always ellipticolij. The in- 
vestigation of a general MHD equilibrium becomes much 
more complicated when one tries to include the effects 
of other important factors, say of viscous stress tensor. 
Recently Ren et al.li3 have studied the deformation of 
magnetic island by including the effect of sheared flow 
and viscosity into an ideal two-dimensional MHD equi- 
librium configuration. However, there is an element|-0|f 
inconsistency, whether an ideal equilibrium is rcalistidI3. 
Heuristically, one ignores the resistivity in the Ohm's law 
while calculating the equilibrium, but then a resistive 
stability analysis based on a stationary equilibrium re- 
mains questionable_as long as the field diffusion is not 
taken into accounto. Montgomery et al.Ej have investi- 
gated the problem on non-ideal static axisymmetric equi- 
librium. There have also been attempts to calculate resis- 



tive axisymmetric equilibrium with only toroidal flowE: 
It has been further argued-that tokamak equilibrium flow 
is cither purely toroidafltl or the poloidal compoiiG|nt is 
smalla and quickly damped by magnetic pumpingcJ. So 
there is a natural tendency to exclude the poloidal flow 
while calculating an equilibrium. But when one consid- 
ers finite conductivity with purely equilibrium toroidal 
flow, the conductivity (hence the resistivity) becomes a 
function, rtf space. In general, the resistivity is not a flux 
functionEHrlij irrespective of equilibrium flow. In this re- 
port, we ask the very pertinent question, whether the sit- 
uation changes in presence of poloidal flow. As we show 
that a uniform resistivity profile is consistent in presence 
of poloidal flow, whereas it has been shown that a scalar 
pressure equilibrium can not have uniform resistivityEj. 
Further we show that in a field-reversed (FRC) axisym- 
metric configuration with no toroidal magnetic field, both 
toroidal and poloidal equilibrium flows can coexist with 
finite resistivity, .which is not found to be the case with 
ideal equilibriumliJ. 

We consider the equilibrium resistive MHD equations 
with plasma flow. The equations are 



V • (pv) = 0, 
V-B = 0,VxE = 0,VxB = 0, 
E + v x B = j/cr, 
p(v • V)v = j x B — Vp, 



(1) 
(2) 
(3) 
(4) 



where the symbols have their usual meanings. We use 
a right handed cylindrical system (r, 6, z) with z as the 
axis of symmetry, 9 as the toroidal angle, and r along the 
major radius of an axisymmetric device. We assume the 
plasma flow to be arbitrary (toroidal and poloidal) and 
axisymmetry is assumed i.e. d/ d6 = 0. The plasma resis- 
tivity 7] = o~ l is assumed to be an unspecified function 
of r and z. We have further assumed here that the equi- 
librium is maintained in a steady-state through resistive 
diffusion. The magnetic induction equation allows us to 
write the magnetic field as 



B = -Vip x eg 

r 



I. 
-eg, 



(5) 
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where ip is the magnetic flux function which is the az- 
imuthal component of the vector potential A and I is 
the current function. Similarly, following the continuity 
equation, Eq.(Q), we can express the plasma equilibrium 
velocity as 



v = — Vp x eg + ujreg, 
pr 



(6) 



where ip is the velocity stream function and uj — vg/r is 
the toroidal angular velocity. We also assume that the 
flow is incompressible i.e. V • v = 0. 

Because the flow is now in both toroidal and poloidal 
direction, the poloidal component of current, j p need not 
vanish. In general the current can be expressed as 



-A*ipeg + -VI x eg, 
r r 



(J) 



where A* is the elliptic operator defined by A*^ = 
r 2 V ■ (^iVip) ■ Taking curl of the Ohm's law Eq.(||), we 
have 



V x (v x B) = V x (j/a) 



(8) 



with 



v x B 



1 



I 



VtpxVip 7 Vip + ujVip. (9) 



pr* pr* 
The eg component of Eq.(0) can be now written as 



e g ■ Vu> x Vtp — eg ■ V — =■ x Vtp 
1 /1„ „ T 2dl „ 2 " 

= — -Vcr - VI + -— - V 2 I 

or \<j r or 



(10) 



We now invoke the large aspect ratio expansion and 
assume that toroidal magnetic field, to the first approx- 
imation can be written Bg ~ B^r^/r. Here, Bq is the 
value of the toroidal magnetic field at center of the cylin- 
drical cross section of the torus at distance ro from the 
axis of symmetry. To this effect we have the current 
function / rs Boro = const. Under these assumption, the 
above equation can be written as 



eg ■ Vuj x Vtp 



■^-^eg -VpxVtp 
p*r* 

21 „ 

ireg ■ Vr x Vtp. 

pr 6 



(11) 



The first term on the right hand side of the above equa- 
tion vanishes by virtue of the incompressibility condition 
and the continuity equation. We neglect the second term 
because of its 1/r 3 dependence and find that the toroidal 
angular velocity vg/r = cu{ip) becomes a surface quantity. 
This also further means that j = jgeg with j p = 0. We, 
however, note that the condition uj = Lj(tp) is identically 
satisfied in a field-reversed configuration (FRC) where 



1 = and large aspect ratio approximation is not re- 
quired. It can be noted here that without any approxima- 
tion, uj becomes a flux function when one considers ideal 
cquilibriumcl or resistive equilibrium with only toroidal 
flowta. 

Now, we consider the momentum equation Eq.(Q) and 
its eg component. With the above approximations, we 
can write Eq.(||) as 



jgeg x Bp = VP + pV 



2p 2 r' 



■(Vtpf 



pr* 



Vtp ) Vtp - Lu'(V<p x Vtp) 



2lu dp 
r oz 



—— eg ~ lu rVr 



(12) 



In the above equation B p is the poloidal component of 
the magnetic field and (') denotes derivative with respect 
to ip. Taking the eg component of the above equation, 
we have, 



e ■ Vp x V(ujr 2 ) = 0, 



(13) 



which means that tp = p(ujr 2 ). We take the simplest 
situation of p oc tor 2 which yields another surface quan- 
tity, tp/r 2 = However, it is important to note that, 
C(^) is not an arbitrary function in the sense that it is 
proportional to the toroidal velocity uj(ip) i.e. the toroidal 
and poloidal flows are no longer independent. Physically, 
one can understand this by noting that finite resistivity 
allows plasma motion across the the flux surfaces. Equiv- 
alently, toroidal flow, in a resistive axisymmetric plasma, 
is always associated with poloidal flow. 

Because of equilibrium flow, however, plasma pressure 
p is no longer a flux function now. Taking the B p com- 
ponent of the momentum equation (Eq.(H3)), we have 



B, 



— + v 



2p 2 r- 



■(Vtpf 



1 2 2 

—tor 
2 



^V-(^V(Cr 2 ))B ; ,VK^. 



(14) 



Depending upon the equation of state, now, several 
options are possible. However, we note that density, in 
general, is not a flux function in presence of arbitrary 
plasma flow, this can be easily seen from the equation 
of continuity Eq. ([!]) , after applying the incompressibility 
condition, 



eg ■ Vp x Vp = 0. 



(15) 



It can be seen from the above expression that p is not a 
surface quantity. We note here that axisymmetric equi- 
libria with incompressible equilibrium flows are genjeijally 
associated with constant density magnetic surfacesli3~EI 
One is also free to choose density as a flux function in 
case of resistive axisymmetric equilibrium with only in- 
compressible toroidal flowcR 
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Taking the eg component of Ohm's law Eq.(g) along 
with Eq.(§), we have an expression for plasma conduc- 
tivity, 



a E r + -r£B, 
P 



A*V = 0, 



(16) 



where Eq is the longitudinal externally applied electric 
field at major radius r — r*o- We immediately see from 
the above expression that conductivity, in general, is a 
space dependent quantity. 

In what follows, we shall consider two cases with 
(?) uniform and constant density and (ii) a nonuniform 
density. In the second case, we consider isentropic mag- 
netic surfaces. We now assume that plasma density is 
uniform and constant i.e. p = const, and normalize our 
equations to p = 1. We can now write Eq.(|l^) as, 



B p V 



1 



A*(Cr 2 )B p -V(Cr 2 ). 



(17) 



Integration of the above equation yields the equivalent 
Bernoulli's equation, 



P + 



2r 2 



+ (18) 



where the integration is along a magnetic field line and 
x(V>) is an arbitrary surface quantity. The solubility con- 
dition further requires that 



-^lA*( C r 2 )B p .V ( Cr 2 ) 



0. 



(19) 



We further assume that part of the pressure gradient 
that varies— .within a magnetic flux tube has no \7ip 
componentEZI i.e. 

V ^ • V j §~ p A * (C^)B P • V(Cr 2 ) = 0. (20) 

Together with Eq.(p"8|) and the above assumption, the Vip 
component of the momentum equation yields the equiv- 
alent Grad-Shafranov equation, 



A*ip + r 2 (x + Juo 2 r) + 



iWl 2 



0, (21) 



with two arbitrary flux functions x(V0 and u(ip)- The 
primes refer derivative with respect to ip. 

We now consider the second case where we consider 
a nonuniform density. With incompressible flow, mag- 
netic surfaces with constant entropy is quite a reasonable 
approximation in ideal MHBd'lII. However, considering 
long resistive diffusion time, the right hand side of Eq.(jJ|) 
can be neglected and we can continue to proceed with 



isentropic magnetic surfacest3. The equation of state can 
now be written as, p = Sp 1 , where, S(ip) is the entropy 
which is a flux function and 7 is the ratio of specific heats. 
We now write B p • Vp/p as B p • V[7Sp 7_1 /(7 — 1)], so 
that equivalent Bernoulli's equation can be written as, 



e(v) 



|--V M V(Cr 2 ))B p .V(Cr 2 ) 
B p p \pr z ' 

7 c>7-l 4- _ 



7- 



1 



2p 2 r- 



-u?r\ (22) 



where Q(ip) is arbitrary. As we have assumed previously, 
it requires a solubility condition and the equivalent to 
the assumption (p0|). We can then continue to write the 
equivalent Grad-Shafranov equation by taking the VV> 
component of the momentum equation and applying the 
Bernoulli's law Eq.(|22|), 



A> 



2 6' + uo'uo 2 r - S' 



,7-1 



— V(Cr 2 
pr z 



7-1 



Wtp ■ V«r 2 ). 



(23) 



In the above equation we have four arbitrary surface 
quantities i.e. Q(tp), co(ip), and S(ip) and the primes de- 
note derivative with respect to ip. 

We have derived the differential equations, equivalent 
to the Grad-Shafranov equation, for resistive axisymmet- 
ric plasma with arbitrary equilibrium flows. These equi- 
librium equations Eqs.(|2l], |23|) have to be solved subject 
to conductivity constraint Eq. (|l6|) . Further, in a field- 
reversed configuration (FRC) with no toroidal magnetic 
field, it can be seen from Eq.(|l^) that both poloidal and 
toroidal flow can coexist. 

We now show that a uniform conductivity profile is 
consistent with resistive axisymmetric equilibria with ar- 
bitrary flow. A simple examination of Eq.(^), though 
reveals that uniform conductivity may be possible with 
scalar pressure equilibrium in presence of flow, it how- 
ever provides no easier way of proving it. We note that 
the usual procedure for solving Eqs.(|2l]) and (^3|) requires 
specifying a priori dependence of the respective arbitrary 
functions on ip. However, in the presence of finite resis- 
tivity, the resistivity constraint Eq.(|l6|) can be used to 
solve for ip, which is uniquely determined if the right 
hand side of Eq.(|l6|) is specifiedt3. It should be noted 
here with caution whether the resultant solution for ip 
corresponds to realistic profiles for other physical quan- 
tities such as pressure, density, velocity etc. However, 
our sole aim, here, is to demonstrate the existence of a 
solution consistent with uniform resistivity in presence of 
flows. 

From Eq.(|l5|) we know that p = p(<p), and assume that 
p oc ip. We now assume that conductivity is uniform in 
space so that the resulting Eq.(|l^) can be written as, 



A*ip- 



— — = (3 

r 2 dz ' 



(24) 
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where a and (3 are arbitrary constants. It is worthwhile 
mentioning at this point that Eq.(p4|) can not be used 
in case of very small resistivity. In the limit of vanishing 
resistivity (large (3 in the above equation), the solution of 
Eq.(j24|) contains short scale spatial dependence (bound- 
ary layers), not present in case of ideal equilibrium and 
may lead to unphysical results. 

Note that the above equation is a elliptic equation and 
can be treated as boundary value problem. Following 
Zheng et air 2 ], we assume a solution of the form 



n=0,l,2,... 



(25) 



for the homogeneous part of Eq. (|24|) . We however re- 
tain the odd terms in the summation to take care of 
the asymmetric-term in Eq.(^4j). For simplicity we as- 
sume that f n (r) = for n > 3, which, however, can 
be extended up to any number of terms if required, 
about which we shall make a comment later. Substi- 



tuting En. (25) in the equivalent homogeneous equation 
for Eq.(J2J) we can solve for the functions Jo. l, 2(f)- The 
homogeneous solution of Eq. ( p4[ ) is then given by, 

ip h (r,z) = air 2 {4r 2 + 16z 2 - a 2 [4(lnr) 2 - 41nr + 2 
-r 2 ] + 8az(r 2 - 2 lnr)} + a 2 [2r 2 (2 lnr 
- 1) - a 2 In r (In r + 1) + Aaz In r + 4z 2 ] 
+ a 3 r 2 [a;(21nr - 1) + 4z] + a 4 r 2 , (26) 

where a^s are arbitrary constants to be determined from 
the boundary conditions. A particular solution of Eq.(|24|) 
is ip r> = /3r 2 (21n?' — l)/4. So the complete solution of 
Eq.© is 



(27) 



which can be verified by direct substitution. For a con- 
ducting circular boundary of a toroidal axisymmetric de- 
vice, the constant flux (ip) contours are shown in Fig.l 
(a) which shows a scaler pressure equilibrium. The solu- 



tion for a 



is shown in Fig.l (b). Note that a oc 



is the only possible solution for resistive axisymmetric 
equilibrium without flowt9. 



In principle the expansion in Eq.(Pq) should be re- 
tained with a large number of terms which will result 
a equally large number of arbitrary constants for the so- 
lution in ip. These constants can then be used to shape 
any arbitrarily shaped plasma boundary. 

In passing, we would like to note that resistive field 
diffusion (dB/dt ^ 0) is intrinsically involved with non- 
stationary equilibria (v ^ 0). However, a series of ideal 
quasi-stationary equilibrium states can be built up with 
dB/dt = in which, the effect of finite resistivity is 
only .to slowly evolve the equilibrium in a diffusive time 
scaleM 
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